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Icing alters the shape and surface characteristics of aircraft components, which results in altered aerodynamic 
forces and moments caused by air flow over those iced components. The typical effects of icing are increased drag, 
reduced stall angle of attack, and reduced maximum lift. In addition to the performance changes, icing can also affect 
control surface effectiveness, hinge moments, and damping. These effects result in altered aircraft stability and 
control and flying qualities. Over the past 80 years, methods have been developed to understand how icing affects 
performance, stability, and control. Emphasis has been on wind-tunnel testing of two-dimensional subscale airfoils 
with various ice shapes to understand their effect on the flowfield and ultimately the aerodynamics. This research has 
led to wind-tunnel testing of subscale complete aircraft models to identify the integrated effects of icing on the aircraft 
system in terms of performance, stability, and control. Data sets of this nature enable pilot-in-the-loop simulations to 
be performed for pilot training or engineering evaluation of system failure impacts or control system design. 
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Nomenclature 

wing span 
drag coefficient 
aileron hinge moment 
aileron control effectiveness 
lift coefficient 

pitching moment coefficient 

pitching moment coefficient derivative with 

AoA, longitudinal stability 

normal force coefficient 

chord length 

ratio of protuberance height to chord length 
Reynolds number 
angle of attack, deg 
angle of sideslip, deg 

control surface deflection: aileron, elevator, 
rudder, flap, deg 


I. Introduction 

A IRFRAME icing has been a threat to aviation safety for nearly 
80 years. Over that time, significant progress has been made to 
reduce the hazard through the design and implementation of ice 
protection systems, improvement of certification procedures, and 
improved pilot training. The foundation of these advancements 
has been research to better understand the icing environment, the 
physics of the ice accretion process, and the degrading aerodynamic 
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effects caused by ice. However, even with these advancements, 
icing-induced loss of control incidents and accidents have occurred 
and continue to occur on all classes of aircraft, from general aviation 
airplanes [1,2] and business jets [3,4] to transport category aircraft 
[5-1 1], These events typically result from a combination of causal 
factors, but they usually precipitate from the altered aerodynamics 
caused by icing. 

The continuing problem of icing has been recognized by many 
organizations including the National Transportation Safety Board 
(NTSB), the Commercial Aviation Safety Team (CAST), the Federal 
Aviation Administration (FAA), and the National Aeronautics and 
Space Administration (NASA). The NTSB has made numerous 
safety recommendations after icing related incidence and accidents 
and has ranked icing among its “MOST WANTED Aviation Trans- 
portation Safety Improvements’’^ since 1997. CAST, a consortium of 
industry and government organizations working to improve aviation 
safety, has developed interventions specific to icing. They have 
also specified that the implicit icing effects on flight envelope 
protection and high-fidelity simulators for upset recovery training are 
among their highest ranked intervention needs [12]. The FAA and the 
NASA Aeronautics Research Mission Directorate have invested 
significant resources to support these safety recommendations, 
develop interventions, and to advance the state of knowledge on the 
icing atmosphere, the ice accretion process, and aerodynamic effects 
of icing. 

The purpose of this paper is to review the methods currently 
available to model and simulate icing effects on performance, 
stability, and control. The primary method to model icing effects use 
wind-tunnel testing of 2-D and 3-D models to develop mathematical 
expressions or databases of aerodynamic forces and moments. Icing 
effects simulation uses these mathematical models or databases in 
conjunction with the equations of motion to explore the inter- 
relationship of various parameters and the modified flight envelope 
imposed by the iced aerodynamics. Icing effects modeling and 
simulation methods have varied maturity levels and limitations in 
terms of verification and validation. Recommendations are made to 
continue research in flight simulation modeling and real-time 
modeling to continue advancing the state of the art. 


s http://www.ntsb.gov/Recs/mostwanted/aviation_issues.htm 
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II. Airframe Icing: How it Forms 

To model and simulate the degraded aerodynamics that results 
from ice, one should understand the ice accretion process and the 
resulting ice characteristics because they drive the aerodynamic 
effects. 

Ice forms on an airplane as it flies through clouds consisting of 
supercooled liquid droplets (i.e., liquid water drops that are below 
0°C). The liquid droplets impinge on the leading edge surfaces, 
exchange heat with the environment and aircraft surface, and change 
phase from liquid to ice. Initially, the ice forms as a thin, rough layer. 
Overtime, the ice thickness increases and reshapes the leading edges 
of all frontal surfaces. The rough, reshaped leading edges, in turn, 
degrade the aerodynamics of the airframe. 

The ice accretion characteristics (size, shape, extent, roughness, 
and translucency) are related to the atmospheric environment param- 
eters, the flight condition, the geometry of the aircraft component 
being iced, and the time in icing conditions. The icing cloud atmo- 
sphere has been defined using three primary parameters. 

1 ) Liquid Water Content (LW C) : a measure of the amount of liquid 
water in a unit volume of space. 

2) Median Volumetric Diameter (MVD): a representative size of 
the water droplet spectrum. MVD is a value for a given spectrum such 
that half the volume of water is in drops with diameters smaller than 
the MVD and half the volume of water is in drops with diameters 
greater than the MVD. 

3) Static Temperature (Ts): the ambient air temperature. 

The flight conditions consist of true airspeed, angle of attack, and 
altitude. Geometric parameters consist of the size of the leading edge 
radius, the chord length (c), single or multi-element wing, and 
straight or swept wing or tail. Each of the preceding parameters 
influences the ice accretion development. The result is a multitude of 
possible ice shapes. 

To illustrate the resulting geometry change due to icing conditions, 
two ice shapes are shown in Fig. 1. They were accreted on a NACA 
23012 airfoil model in the NASA Icing Research Tunnel [13]. 
Figure l_a shows an example of a horn ice shape that forms at tem- 
peratures near freezing, where impinging water droplets can flow 
before freezing. Figure fb shows an example of a streamwise ice 
shape. This type of ice forms at temperatures well below freezing 
where all of the impinging droplets freeze immediately on contact 




Fig. 1 Ice shapes accreted at different conditions in the NASA Icing 
Research Tunnel: a) Ice accreted at temperature near freezing, b) Ice 
accreted at temperature well below freezing. 


with the airfoil. Because these shapes are so different, they have very 
different effects on the aircraft aerodynamics. 

To reduce the aerodynamic penalty of icing, ice protection systems 
(IPS) [14] are typically incorporated into airplanes to prevent or 
remove ice from critical flight components such as the wing and tail, 
engine nacelles and air intakes, propellers, and others. Ice protection 
systems that prevent ice from forming are typically called anti-icing 
systems, whereas those that remove ice are called de-icing systems. 
Anti-icing systems typically use heat from engine bleed-air or 
electrothermal heaters to evaporate water drops that impinge on the 
leading edge surfaces or to prevent the water drops from freezing on 
the heated areas. De-icing systems debond ice that has already 
formed on the leading edge through a mechanical deformation of the 
surface (e.g., pneumatic de-icing boots) or through intermittently 
heating the surface (e.g., electrothermal heaters near the propeller 
root). The design of an IPS on any given airplane may use a com- 
bination of anti-icing and de-icing equipment. Some airplanes have 
bleed-air thermal anti-icing on the wings and engine nacelle leading 
edges only, leaving the tail surface completely unprotected. Others 
have pneumatic de-icing boots on the wing, and on the horizontal and 
vertical tails while having electrothermal de-icing on propellers. 
Some have a combination of bleed-air thermal anti -icing on the wing 
and engine nacelle leading edges, pneumatic de-icing boots on the 
horizontal tail, and no-ice protection on the vertical tail. 

There are other IPS designs as well, so the list could go on. In the 
end, the airframe manufacture determines the locations in which ice 
protection is needed and what type of IPS is required to enable 
safe flight operations in icing conditions and meet certification 
requirements. 

Airframe icing is a highly complex dynamic phenomenon. The ice 
that forms on an airplane is the result of a combination of icing cloud 
atmospheric parameters, flight conditions, airframe geometry, loca- 
tion and type of ice-protection systems, and the time flown in icing 
conditions. Icing atmospheric conditions are not static but change 
spatially and temporally. Similarly, flight conditions and airplane 
geometry change as the phase of flight changes from climb, cruise, 
descent, to approach and landing. Some ice-protection systems do 
not remove all ice, but leave residual ice on the protected surface. 
Some thermal IPS may not fully evaporate the impinging water and 
allow water to run back aft of the protected areas where it freezes to 
form frozen rivulets, an ice ridge, or other shapes. The amount of time 
in icing conditions also plays a significant role as ice builds up on the 
unprotected surfaces such as unprotected tail, struts, radome, engine 
pylons, and flap hinge fairings. 

Clearly, the inputs to icing formation are multidimensional and 
result in a wide range of potential ice shapes and ice characteristics. 
This presents a large challenge to accurately predict the entire range 
of possible effects on aircraft performance, stability and control, and 
handling. Over the past 80 years, many efforts have been undertaken 
to advance the understanding of ice accretion and the impact on 
aerodynamics. The next section briefly discusses some of these 
advances. 


III. Modeling Icing Effects: Wind-Tunnel Results 
From Two-Dimensional Wing Models 

The typical aerodynamic effects of ice on airfoils are increased 
drag, a reduced stall angle of attack, and a reduced maximum lift. In 
addition to these primary parameters, icing affects the pitching 
moment and, if equipped with a control surface, the control effec- 
tiveness and hinge moment. 

To understand and model these effects, research organizations 
have conducted numerous wind-tunnel tests on a number of 2-D 
airfoil models and some 3-D wing and tail models. Recent iced 
aerodynamics simulation studies provide a significant body of 
knowledge examining the effects of ice accretion (size, shape, 
roughness, 3-Dity) on drag increase, lift loss, and stall characteristics 
of a NACA 23012 airfoil at full-scale and subscale geometry and 
Reynolds numbers [13]. Figure 2 illustrates the change in perform- 
ance characteristics with different ice shapes at Re = 1.8 x 10 6 . 
Increase in drag is observed at the initial onset of icing (roughness), 
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Fig. 2 Aerodynamic performance comparison for NACA 23012 model with casting simulations (Ref. [13]). 


and drag continues to increase as the ice takes shape and continues to 
grow in size. Drag increase is observable at all angles of attack. Ice 
affects lift mostly at higher angles of attack in which early stall 
occurs. For roughness, this is due to early flow separation caused by 
loss of boundary -layer momentum. For larger ice shapes, early stall 
occurs because the separation bubble that forms aft of the ice shape 
cannot reattach to the airfoil. Before this new stall point, there are 
shifts in the lift curve slope and the stall break, and post-stall behavior 
is clearly affected by the ice. The pitching moment is also affected 
due to the altered pressure distribution over the airfoil, and a pitch up 
tendency occurs at significantly lower angles of attack compared 
with the clean (no-ice) baseline. 

Ice accretions have also degraded control surface effectiveness and 
increased the hinge moment loads on 2-D airfoil models. Figure 3 
shows the effect of a simulated ice shape on aileron effectiveness 
(C L g A ) and hinge moment (C A M ). This data are for a forward facing 
quarter round, k/c = 0.0139 at various chordwise locations at Re = 
1.8 x 10 6 [15], The first plot shows the rate of change in lift 
coefficient with change in aileron deflection, and the second plot 
shows the rate of change in hinge moment with change in aileron 
deflection. On the clean airfoil, the C ! M value remained relatively 
constant from AoA = —5 deg to 5 deg. As the angle of attack was 
increased from this point, there was a gradual, almost linear, reduc- 
tion in aileron effectiveness to AoA = 14 deg where the airfoil 
stalled. At stall, there was a sudden reduction in aileron effectiveness, 
where C L $ A became negative, meaning that deflecting the aileron 
trailing edge down at this angle of attack decreased lift. All three ice- 
shape locations (xj c = 0.02, 0. 10, and 0.20) showed reduced aileron 
effectiveness at AoA = 0 deg, when compared with the clean air- 
foil. However, the reductions were greater when the ice shape was 
located at x/c = 0.10 and 0.20 than when it was located at 
x/c = 0.02. For all three iced cases, the aileron effectiveness started 
to decrease at lower angle of attack than the clean model. However, 
the rate at which it decreased varied with ice-shape location. The 
hinge moment data did not show any noticeable trend in the rate of 
change in hinge moment with aileron deflection. The C h SA values did 
not appear to vary significantly with change in angle of attack. Also, 
the presence of ice shapes did not appear to have altered the C h $ A 
values significantly either. 


Experimental efforts such as these have enabled a comprehensive 
understanding of the basic flow physics associated with various ice 
accretions and helped identify the major aerodynamic penalties 
associated with the particular ice accretion studied. Bragg et al. [16] 
provided an in-depth analysis of flowfields caused by ice geometries 
categorized as roughness, horn ice, streamwise ice, and spanwise 
ridge ice. Studies such as these are critical in developing and 
validating computational fluid dynamics (CFD) tools to predict iced 
airfoil and eventually iced airplane, aerodynamic characteristics. 
These studies are also beneficial in defining ice-shapes character- 
istics that are critical to specific aerodynamic parameters. However, 
to understand the interaction of iced wings and tails for various ice 
shapes on airplane flight dynamics, a different approach and test 
methods need to be employed. 

IV. Modeling Icing Effects: Wind-Tunnel Results 
From Subscale, Complete Airplane Models 

When considering the aerodynamic effects of ice on an airplane 
configuration, the effects found on 2-D airfoils apply, but are 
expanded upon. In this case, the effects of drag on the airfoils are not 
isolated, but couple into the pitching moment of the airplane. The lift 
degradation on the horizontal tail also couples into the airplane 
pitching moment in a static sense but also a dynamic sense by reduc- 
ing pitch damping. Lift degradation on the wing can affect longi- 
tudinal stability and roll damping. Adding the consequence of 
reduced control surface effectiveness can lead to an aircraft with 
substantially degraded stability and control and handling qualities. 

These coupling effects are dependent on the airplane design, and 
airplane manufacturers are required to demonstrate safe flying 
characteristics with potential ice accretions that may occur in both 
normal operation and failure of the ice-protection system. To manage 
the wide range of possible ice shapes as described in Sec. II, the ice 
shapes tested are defined based on flight scenarios, but the onus is on 
the manufacturer to define the most critical ice shape for the given 
scenario. The work that has been and continues to be performed in 
aerodynamic studies such as [1_6] helps identify the critical ice-shape 
features and verify applicability of data from subscale testing. In 
practice, the ice shapes now commonly tested are 1) ice roughness 
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Fig. 3 Aileron effectiveness and hinge moment on NACA 23012 with simulated ice shapes at various chordwise locations (Ref. [15]). 
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(simulating delayed IPS activation or residual/intercycle ice), 
2) runback ice (ice that forms aft of thermal IPS), 3) failure ice (ice 
that forms on the IPS after a failure, typically large horn ice), and 
4) holding ice (ice that forms on unprotected surfaces during a 45 min 
hold, typically large horn ice). 

Two test methods have been used to collect data on subscale, 
complete airplane models with ice shapes. These are static testing 
and static and dynamic testing. Each of these will be discussed in the 
following subsections. 

A. Static Testing of Subscale, Complete Airplane Models 

Static testing produces databases of airplane force and moment 
coefficients for static conditions (a, /3, 5 control surface , ice case). These 
databases provide valuable insights into the effects that ice has on the 
airplane stability and control. Two examples of this type of modeling 
are provided to demonstrate this. 

Prompted by icing related incidents, C. L. ’Kelly’ Johnson 
investigated the effects of ice on performance, stability, and control 
on the Lockheed Electra [17]. He found through wind-tunnel tests 
using artificial ice shapes with a scaled model (Fig. 4) that maximum 
lift was reduced by 32%, drag was increased by 47%, stall angle of 
attack was reduced by 4°, and aileron control effectiveness was 
reduced by 36%. Figure 5 shows the effect of icing on roll control. 
The results were obtained at Re = 1.6 x 10 6 . These data clearly 
show the dramatic reduction in stability and control when the Electra 
was iced up. 

More recently, NASA conducted a static wind-tunnel test on a 
twin-engine short-haul transport to measure icing effects on perfor- 
mance and stability [18] at the NASA Langley Research Center 
(LaRC) 14 by 22 ft wind tunnel. These tests were conducted on a 1 78- 
scale model (Fig. 6) at Re = 1.7 x 10 6 . Representative results from 



Fig. 4 Subscale Lockheed Electra with artificial ice shapes (Ref. [171). 



Fig. 5 Effect of ice on aileron control: full aileron deflection (Ref. [171). 



Fig. 6 NASA Langley Research Center 1 / 8-scale twin-engine subsonic 
transport model. 


this test are shown in Fig. 7. In this case, the ice shapes on the wings 
reduced stall angle of attack by about 5 deg and lowered maximum 
lift coefficient (C Lmax ) by approximately 16% compared with the 
clean data at the test Reynolds number. The difference between the 
iced configuration and the no-ice configuration at flight Reynolds 
number is expected to be greater due to increase in C L max for the no- 
ice configuration. A positive shift in the pitching moment occurred 
with the ice cases, but longitudinal stability ( C ma ) was approximately 
the same as the no-ice case until the AoA reached about 8 deg where 
stall occurred and the aircraft became longitudinally unstable. No 
data were acquired with control surface deflections other than flaps in 
this wind-tunnel test, so that control effectiveness with the ice shapes 
could not be evaluated. 

B. Static and Dynamic Testing of Subscale, Complete Airplane 
Models 

Although static testing provides insights into the effect ice has on 
stability and control, it does not address the effects on dynamic 
motion. To better understand and model the full range of an iced 
aircraft flight dynamics, data from forced oscillation and rotary 
balance tests need to be gathered. 

To that end, static and dynamic wind-tunnel tests were conducted 
under various projects on scaled DeHavilland DHC-6 Twin Otter, 
Cessna business jet [19], and Lockheed S-3B Viking models at the 
Bihrle Applied Research Large Amplitude Multi-Purpose (LAMP) 
facility in Neuburg, Germany (Fig. 8). The scale of the models were 
6.5, 8.3, and 6.5% and they were tested at Re = 0.13, 0.15, and 

0. 20 x 10 6 , respectively. 

The primary purpose of these research efforts was to understand 
the effect of airframe icing on flight dynamics. This was accomp- 
lished through the creation of flight simulation models that used the 
aerodynamic databases derived from wind-tunnel test results. This 
process was used to explore the utility of iced flight simulation in 
pilot training applications using the Twin Otter aircraft. Simulation 
models were developed for the no-ice (clean) and two IPS failure 
ice-shape configurations. The successful implementation of this 
approach led to a similar effort using a Cessna business jet. In this 
case, four configurations were considered: no-ice (clean), ice rough- 
ness, wing IPS failure ice, and wing runback ice. Lastly, since NASA 
John H. Glenn Research Center at Lewis Field is modifying a 
Lockheed S-3B Viking for icing flight research, there is a need to 
understand the impact of potential ice accretions on its flight 
dynamics. In this case, data were acquired for the no-ice baseline 
configuration and icing configurations consisting of large horn ice 
shapes on the wing and tail leading edges to represent a hold in icing 
conditions, and spanwise ridge shapes aft of the thermal IPS on the 
wing and horizontal tail to represent runback ice. 

1. Challenges Associated With Subscale, Complete Airplane Model 
Testing 

The geometric scale of these models and the associated Reynolds 
numbers required the development of methods to account for 1 ) the 
premature stalling characteristics typical of low Re and 2) ice 
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Fig. 7 Effects of ice on longitudinal aerodynamics of twin-engine transport Sr = 40 deg (Ref. [18]). 


accretion scaling to represent full-scale iced aerodynamics. To 
develop flight simulation models that are representative of full-scale 
flight characteristics from the low -Re wind-tunnel data, the clean no- 
ice wind-tunnel data are projected out along a trend line to the angle 
of attack in which stall is anticipated at flight Reynolds numbers. This 
angle of attack can be determined from other sources such as high 
Reynolds number wind-tunnel data and flight tests. This method of 
data extension is shown in Fig. 9 for the Twin Otter and Fig. 10 for the 
Cessna business jet. Another challenge with this type of wind-tunnel 
data is the reduced effectiveness of flaps for high lift (Fig. 9). As flap 
angle increased, the change in lift associated with that flap deflection 
was not as great as seen in flight. In this case, corrections were made 
based on flight data. 

The small geometric scale of the complete airplane test articles 
required methods for defining the size and position of artificial ice 
shapes to represent full-scale aerodynamics. Although research on 
subscale 2-D airfoils with large ice shapes has shown iced aerody- 
namics to be relatively insensitive to Re effects [16], these complete 
aircraft model tests were conducted at much lower Re, and further 
research was warranted. To that end, a series of wind-tunnel tests on 
full-scale, midscale, and small-scale test articles were conducted. 
The small-scale test articles were representative of the subscale 
complete aircraft models. For the Twin Otter, this consisted of testing 
a full-scale, a 42%-scale and a 7%-scale horizontal tail (Fig. IT) to 


identify the size and position of the artificial ice for the 6.5%-scale 
Twin Otter model [20]. The full-scale ice shape for the wing and 
horizontal and vertical tail are shown in Fig. 1_2. 

For the Cessna business jet, these tests used a full-scale, a 41. 7%- 
scale and an 8.3%-scale wing panel model (Fig. 13) to identify the 
artificial ice for the 8. 3%-scale business jet model [21_] . The full-scale 
ice shapes that were tested are shown in Fig. 1_4. 

These tests provided valuable insight into the aerodynamic scaling 
relationships for artificial ice shapes. For example, the failure IPS ice 
shapes could be geometrically scaled to obtain representative aero- 
dynamics with the subscale model at low Re conditions. The 
roughness and runback ice cases for the Cessna business jet proved to 
be greater challenges. For these ice shapes, geometric simulation did 
not produce acceptable results because the full-scale C Lmax with ice 
shapes were greater than the subscale no-ice baseline C L max . This led 
to determining the size and position of artificial ice on the subscale 
model that had similar offsets to those observed in the full-scale tests. 
Those results were then projected to the flight Re condition similar to 
the no-ice baseline simulation models. 


2. Aerodynamic Model Results 

The aerodynamic scaling tests provided confidence that the 
measurements made in the LAMP facility could be used for 



Fig. 8 Scale models of DHC-6 Twin Otter, Cessna business jet, and S-3B Viking. 
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Fig. 9 Comparison of Twin Otter no-ice wind-tunnel normal force 
coefficient data with simulation model values for flight Re. 


modeling full-scale flight dynamics. The data acquired at the LAMP 
consisted of forces and moments along all three axes for static 
conditions (a, /J, 8 a . 8 e . 8 r , 8f)\ forced oscillations in roll, pitch, and 
yaw axes; and steady rotation about the wind vector. These data were 
acquired for the no-ice baseline and for the iced configurations. The 
no-ice baseline data were shifted to a full-scale Re stall AoA as 
described earlier. To illustrate key effects that ice had on the 
aerodynamic characteristics, select aerodynamic model results from 
the Twin Otter effort [22,23] are presented. 


C. Normal Force 

The normal force coefficients were measured over a wide range of 
angles of attack, sideslip angles, elevator deflections, rotational 
velocities, and pitch rates during forced oscillations. Data from these 
measurements were analyzed to determine the model structural 
dependencies and to insure the preservation of all nonlinear effects. 
The initial definition of the basic force model ( C\' basic ) was built from 



a, deg 



Fig. 10 Comparison of business jet no-ice C L and C,„ wind-tunnel data 
with simulation model values for flight Re. 
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Fig. 11 Depiction of Twin Otter horizontal tail at various scales. 


the most significant functional dependencies. For the Twin Otter 
these were angle of attack, angle of sideslip, and flap deflection. 
Incremental coefficient data tables were then generated to provide the 
effects of control deflections, dynamic damping, and power effects. 
The normal force coefficient can be represented by Eq. (1_). The 
C;v BASIC lookup table is based on static inputs of a, fi, Sf. The A C N 
lookup table of increments of C N is based on inputs of a, 
Q.b/2V x SGN(/3), \f}\, Sf. The A C Nde lookup table of increments of 
C N is based on static inputs a, 8 e ,8f. Lastly, the A C N lookup table of 
increments of C N is based on inputs a, qc/2V, 8f. The result from 
each of these lookup tables is added to arrive at the total normal force 
coefficient. 

1. Nonna! Force Coefficient Model for No-Ice Baseline 

C A TO talclean = CaTbasic^’ P ' + AC V Rtr r ( a ’ ^ b / 2 V 

x SGN(/3), \P\,8f) + ACjv de (o!, 8 e , 8 f ) 

+ A C Na (fx,qc/2V,8 f ) (1) 

where SGN(/3) = ±1 depending on the sign (±) of f J >. The 
differences in the no-ice and ICE02 (IPS failure ice on wing, 
horizontal and vertical tail) data sets drove the modeling effort to 
create separate databases for each configuration. Equation (2) shows 
that the model structure for ICE02 is the same as the no ice, but the 
values within the tables are different. 

2. Normal Force Coefficient Model for All-Iced ( ICE02 ) 

C a TOTA liceo 2 = i2 + AC a ROTI2 (“- & b / 2 V 

X SGN(j8), \P\,8f) + AC Nmn (a, 8 e , 8 f ) 

+ A C NQ Ja.qc/2V,8 f ) (2) 

As an example of the C, Vbasic lookup tables, the normal force 
coefficient data for the Twin Otter with flaps at 0, 20, and 40 deg are 
presented in Fig. 15 for the no-ice baseline and the ICE02 
configurations. Comparing these results, the effect of ice on the 
normal force is observed mostly at the angle of attack near stall, 
where there is a reduction in the maximum normal force for all flap 
settings. This effect is most pronounced with the flaps setting of 
8f = 20. The reduced normal force is similar to findings with 2-D 
airfoil research, but note that the stall break and post-stall character is 
similar to the no-ice baseline. This small change is due to the rather 
benign stall characteristics of the Twin Otter wing and the size and 
shape of the wing ice tested. 

D. Pitching Moment 

As with the normal force, the pitching moment coefficient data 
were measured over a wide range of angles of attack, sideslip angles, 
elevator deflections, rotational velocities, and pitch rates during 
forced oscillations. Data from these measurements were tabulated 
into five databases so that the pitching moment coefficient can be 
represented by Eq. (3). Similar to the C N equations, there is a 
corollary equation and set of lookup tables for the ICE02 
configuration. There is an additional term in the C m equation for the 
thrust effects: A C m . This database was developed using existing 
NASA Twin Otter flight data. 
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Fig. 12 Full-scale ice-shape profiles for Twin Otter wing, horizontal tail and vertical tail. 


1. Pitching Moment Coefficient Model for No-Ice Baseline 

^'"TOTAL CLEAN = ^'"BASIC P' ^/) ^ ^'"ROT ( “ ’ &b/2V 

x SGN(/5), \P\,8 f ) + AC %B (a, 8 e , 8 f ) 

+ A C mg (a, qc/2V, 8 f ) + A C„, CT (a, C T , 8 f ) (3) 

As an example of the C m databases, the pitching moment 
coefficient data for the Twin Otter with flaps at 0, 20, and 40 deg are 
presented in Fig. 16 for the no-ice baseline and the ICE02 
configurations. The effect that ice has is influenced by the flap 
deflection 8 f . With 8 f = 0 deg, there is a general positive shift in 
C,„ for AoA > —7 deg, indicating the download from the hori- 
zontal tail is reduced for the given elevator deflection. Also, the 
static longitudinal stability ( C ma ) is reduced significantly for AoA< 
—7 deg, indicating longitudinal instability at this range. With 8 f = 
20 and 40 deg, the differences caused by the ice are insignificant in 
the linear range, but reduced static longitudinal stability occurs for 
AoA > 8 deg and for AoA < —4 deg. The wing ice causes the 
reduced C ma at positive AoA, whereas the ice on the horizontal tail 
reduced the C ma at the negative AoA. 

To illustrate the data in the AC,„ DE lookup tables. Fig. F7 is 
provided. Overall, the elevator is effective in creating proper pitching 


moments throughout the AoA range with or without ice. Flowever, 
the ice does reduce the amount of pitching moment affected by the 
elevator in the — 10 < AoA <10 deg range. This means that larger 
elevator deflections are required to trim the iced aircraft. One thing to 
note is this effect is more pronounced with greater flap deflections. As 
flaps are deflected, the inflow angles at the horizontal tail increase 
causing separation bubbles to extend and reducing the effectiveness 
of the elevator. As seen in Fig. 1_8, the effect of the ice is greater in the 
8f = 40 case and extends to lower angles of attack than when the 
flaps are not deflected. 

E. Model Equations for Remaining Force and Moments 

For completeness, the remaining force and moment equation 
models for the no-ice baseline Twin Otter are shown below in 
Eqs. (4-7). Like the normal force and pitching moments, each of 
these equations has identical model structure for the iced case but 
different values within the lookup tables. 

1. Axial Force Coefficient Model for No-Ice Baseline 

c a totalclean = C ABAS JaJ,8 f ) + A C At '(a,8 e 8 f ) 

+ A C Aq (a, qc/2V, 8 f ) + A C AC >, C T , 8 f ) (4) 



Fig. 13 Full-scale, 41.7%-scale, and 8.3%-scale business jet wing panel models. 
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Fig. 14 Ice shapes tested on full-scale business jet wing panel. 


2. Side Force Coefficient Model for No-Ice Baseline 

^V'lOTAI CLEAN ^^RASIC (a,/U f ) + AC YRm (aMb/2V 

x SGNG8), \fi\ ,8 f ) x SGN(/J) + AC r J«, |U<5/) 
x SGN(<5 fl ) + A C 7|r (a, |3 r |) x SGN(<5 r ) 

+ ACy p (a, pb/2V, S f ) + ACy r (o;, /T?/2V, 8 f ) (5) 

3. Rolling Moment Coefficient Model for No-Ice Baseline 

C/total clean = C 'basic(“’ P ' S /) + AC 'rot(“' & b / 2V 

X SGN(0), |/J|, ^ X SGN(/5) + AC (i >, S fl , « 7 ) 

+ AC, ir (a, « r , 5/) + A C lp (a, pb/ 2V, 8 f ) 

+ AC lf (a,ri/2V,« / ) (6) 

4. Yawing Moment Coefficient Model for No-Ice Baseline 

^ n TOTAL CLEAN BASIC p, S f ) + A C„ ROT (a, £2b/2V 

x SGNG6), \ji\,S f ) x SGN(j8) + AC„Ja, |U */) 
x SGN(i5 a ) + A C n Ja, \8 r \,8 f ) x SGN(<5 r ) 

+ AC, lp (a, pb/ 2V, S f ) + A C„ r (a, rb/2V, 8 f ) (7) 

V. Simulating Icing Effects 

A. Implementing Flight Simulation Models 

For the purposes of the research, to understand the icing effects on 
flight dynamics, the flight models described earlier were implem- 
ented using D-Six, (Bihrle Applied Research, Inc.) a commercial 
off-the-shelf product from Bihrle Applied Research, Inc. D-Six is a 


PC-based simulation environment that fully supported both the Twin 
Otter and Cessna business jet simulation models development, as 
well as the analysis and validation activities. 

One key objective of these research efforts was to enable real-time, 
pilot-in-the-loop simulations to demonstrate icing effects on flight 
dynamics to pilots and engineers. D-Six provided the simulation 
environment to accomplish this because it permits the dynamic 
linking of other object modules that can control everything from the 
simulation integration scheme to external graphics and network 
communications with no requirement to edit the source code. Pilot- 
in-the-loop simulations will be further discussed in the Using Flight 
Simulation Models V.Csubsection. 


B. Validating Flight Simulation Models 

New flight test data were required to validate the flight simulation 
models. For the Twin Otter, a flight test program was conducted with 
the no-ice baseline, artificial ice on the horizontal tail only (ICE01 
configuration), and artificial ice on the wing, horizontal tail and verti- 
cal tail (ICE02 configuration) [24]. Flight test maneuvers included 
control doublets, idle -power stalls, steady-heading sideslips, thrust 
transitions, throttle sweep, wind up turns, and simulated approach 
and missed approach. 

Flight data from these maneuvers were then imported into D-Six 
for analysis and validation of the simulation models. A validation tool 
called Overdrive enabled the validation of the simulation aerody- 
namic database against flight-extracted data using the process illus- 
trated in Fig. 19. At each time slice. Overdrive extracts aerodynamic 
moment coefficients from the flight-recorded time histories, as 
shown on the right side of Fig. 19. Angular rates are numerically 
differentiated to obtain the angular acceleration of the vehicle. After 
the removal of the inertial effects, the remainder is nondimension- 
alized to generate the aerodynamic force and moment coefficients 
experienced during flight. 



Fig. 15 Twin Otter normal force coefficient for various flap settings 
with no-ice baseline and ICE02 configurations. 



Fig. 16 Twin Otter pitching moment coefficient for various flap 
settings with no-ice baseline and ICE02 configurations. 
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Baseline, Ice02, 
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Fig. 17 Twin Otter pitching moment coefficient increment for no-ice 
baseline and ICE02 due to elevator deflection (S f = 0). 


Baseline, Ice02, 

§e S e 



Fig. 18 Twin Otter pitching moment coefficient increment for no-ice 
baseline and ICE02 due to elevator deflection (S f = 40). 


At each time step, flight-recorded states, such as angle of attack, 
angle of sideslip, control surface positions, etc., are used to exercise 
the aerodynamic model in accordance with the aerodynamic model 
specification discussed previously. Each aerodynamic model 


element (i.e., pitching moment due to elevator, etc.) is stored and 
summed as prescribed in the aerodynamic model. By over-plotting 
the model predicted coefficients with the flight-extracted total 
coefficients (Fig. 20), differences can be easily identified. Correlat- 
ing the discrepancies with the excitation of the individual elements 
and parameters from the flight time history aids in isolating 
potential weaknesses in the aerodynamic model. 

Overdrive was used with the Twin Otter flight data to validate the 
models. The validation effort was reported on [25] but is currently not 
published. 


C. Using Flight Simulation Models 

One key objective of these research efforts was to enable real-time, 
pilot-in-the-loop simulations to demonstrate icing effects on flight 
dynamics to pilots and engineers. Although the simulation models 
could be run on a laptop or desktop PC, a portable flight training 
device called the Ice Contamination Effects Flight Training Device 
(ICEFTD) (Fig. 21_) was developed to be a more effective tool [26], 

The ICEFTD consists of a raised platform and framework that 
supports a pilot seat, a control yoke, rudder pedals, a twin turbo-prop 
throttle quadrant, a flap selector, three flat-panel monitors for out- 
the-window graphics, and two additional flat-panel monitors for 
instrument panel graphics. The control column is connected to a 
programmable loader for longitudinal force feedback, whereas the 
yoke (lateral) and rudder pedals force gradients are provided by 
spring resistance. Two PCs using D-Six host the simulation models 
and control the graphics. A third PC, the control loading computer 
(CLC), controls the electromechanical loader to simulate represen- 
tative column forces. These PCs are mounted under the floor of the 
device, and the control loader device is mounted to the frame forward 
of the pilot’s feet. A curtain surrounds the ICEFTD to isolate the pilot 
from external visual distractions. The design is well suited for 
mobility and usability at various settings, from labs or offices to class 
rooms or hangars. 

Both the Twin Otter and Cessna business jet flight models have 
been implemented on the ICEFTD and were used for pilot education 
and training, as well as pilot evaluations of the simulation models. 
The ICEFTD has been used to demonstrate icing effects on Twin 
Otter flight dynamics to over 150 pilots at seminars and short courses 
held by the University of Tennessee Space Institute, National Test 
Pilot School, Flight Safety International, and several conference 
exhibits [27]. The ICEFTD was also used in conjunction with the 


a, (3, speed, altitude Linear 

angular rates (P, Q, R) Angular rates accelerations 



Fig. 19 Overdrive process diagram. 


210 


RATVASKY, BARNHART, AND LEE 


2.0 
1.5 
O 1-0 
0.5 
0.0 






Fig. 20 Sample Overdrive result from Twin Otter wing stall maneuver = 0, all iced. 



flight model validation flight tests of the Cessna business jet [28], 
Simulator sessions were conducted before the flight tests to identify 
anomalous flight characteristics that were predicted by the simulator. 
After the flight tests and flight model updates based on the new flight 
records, Cessna flight test pilots reevaluated the simulation models 
on the ICEFTD. For both airplane types, pilots could readily see the 
changes in performance, stall characteristics, and the increased 
workload to recover from the stall with the iced airplane. 

VI. Future Directions 

As described earlier, there is a large assortment of ice shapes 
possible for any aircraft, each having a range of effects on the 
aerodynamics of individual wings, fuselage, empennage, and other 
surfaces. More important, these ice shapes have a wide range of 
effects on the aerodynamics of the aircraft as an integrated machine- 
human system. The research with iced 2-D airfoils has helped reduce 
the scope of ice shapes that need to be considered by identifying key 
ice-shape features and the resulting flowfields. Similar research 



Fig. 21 ICEFTD used for pilot education and training. 


efforts are needed with 3-D wings and tails with control surfaces to 
determine if 3-D ice features such as scallops need to be considered 
as part of the so-called critical ice shapes. 

Regarding icing effects on flight dynamics, the research efforts 
with the Twin Otter and business jet simulation models are an 
excellent start in that they have provided a method to develop and 
validate icing effects simulation models. These flight simulation 
models and the flight training device are valuable resources for addi- 
tional research into icing effects on flight dynamics, pilot education, 
and future flight training simulator requirements. But the knowledge 
gained through the efforts with the Twin Otter and business jet cannot 
be applied across the entire spectrum of aircraft. Further efforts are 
needed to identify icing effects on other classes of aircraft, such as 
regional jets, large transports, and future designs such as blended 
wing body. 

Understanding these needs, NASA is currently using the Twin 
Otter flight simulation models to develop and test new methods for 
identifying icing effects in real time to provide envelope protection to 
avoid loss of control. This effort is being conducted through a NASA 
Research Announcement (NRA) with the University of Tennessee 
Space Institute and Bihrle Applied Research, Inc. 

NASA is also initiating a new icing effects simulation activity 
using the NASA Fangley Research Center Generic Transport Model 
(GTM). A significant database has already been developed using a 
generic twin-engine large transport model configuration to develop 
recovery strategies from loss of control events. To date, the emphasis 
of the effort has been on damaged aircraft or failure conditions, such 
as a rudder hardover. The new icing effects effort will provide another 
scenario under which loss of control can occur. Flight simulation 
models will be developed from this database and used in conjunction 
with adaptive control methods development and testing. 


VII. Conclusions 

The general effects of icing on aerodynamics are increased drag, 
reduced stall angle of attack and maximum lift, and altered pressure 
distribution over the airfoil surface. In addition, the stall break and 
post-stall characteristics can be dramatically different in the iced 
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cases when compared with the no-ice baseline. These aerodynamic 
effects have been studied using 2-D airfoil sections and much has 
been learned regarding ice-shape features and the associated 
flowfields that result in these performance degradations. 

To understand the integrated effects of icing on aircraft perform- 
ance, stability, and control research has been conducted using 
subscale complete aircraft models. From this research, flight simul- 
ation models were developed that incorporate the nonlinear nature of 
icing effects on the forces and moments along all three axes. These 
flight models can be implemented into piloted flight simulators for 
pilot evaluation and loss of control recovery strategies. 

Significant progress has been made in this area of icing flight 
dynamics research, but the knowledge gained with the limited 
number of models may not be applicable to other airplane config- 
urations. Additional research efforts using other airplane models are 
needed to extend our current understanding of icing flight dynamics. 
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